Operation of universal gates in a solid state quantum computer based on clean 
Josephson junctions between d-wave superconductors 
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The operation of solid state superconducting quantum computer based on clean Josephson junc- 
tions between two d-wave superconductors is considered. We show that freezing of passive qubits can 
be achieved using a dynamic global refocusing technique. Further, we demonstrate that a universal 
set of gates can be realized on this system, thereby proving its universality. 

Quantum computation algorithms promise of enor- 
mous speed up in dealing with certain classes of problems 
jl],D can only be realized if a quantum computing device 
is built on a scale of at least several thousand qubits. 
The inherent scalability of solid state devices and high 
level of expertise existing in industrial electronics and 
experimental mesoscopic physics make solid state-based 
quantum computers an attractive choice |SQ . The prob- 
lem of quantum coherence preservation in such devices, 
in the presence of macroscopic number of degrees of free- 
dom, is difficult but at least theoretically solvable 
Moreover, in a recent experiment on a superconducting 
quantum dot (single electron transistor, SET) || coher- 
ent quantum beats were demonstrated in this mesoscopic 
system, which proves its suitability as a qubit prototype. 

The coherent ground state and gapped excitation spec- 
trum in superconductors make coherence preservation 
more achievable; there exist already several suggestions 
for quantum computers based on Josephson junctions 
and superconducting SETs 

In this paper we consider operation of quantum gates 
in a solid state quantum computer based on clean Joseph- 
son junctions between d-wave superconductors (i.e. bal- 
listic DND or D-(grain boundary)-D junctions) Ter- 
minal B of the junction (Fig.l) is formed by a massive 
d-wave superconductor; in a multiple-qubit system, B 
would be a common "bus" bar. Terminal A is small 
enough to allow, when isolated, quantum phase fluctu- 
ations. It is essentially the sign of the superconducting 
phase difference ip between the terminals A and B that 
plays the role of "spin variable" of quantum computing. 
The collapse of the wave function is achieved by con- 
necting terminal A to the equilibrium electron reservoir 
("ground") through a parity key (superconducting SET), 
thus blocking phase fluctuations due to phase-number 
uncertainty relation S. Other parity keys, with differ- 
ent parameters, are used to link adjacent qubits, allow- 
ing for controllable entanglement. (A parity key only 
passes Cooper pairs, and only at a certain gate voltage 
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FIG. I. a) Superconducting DND qubits: A,B are d-wave 
superconductors, N normal conductor, PK parity key, M scan- 
ning tip, il the mismatch angle between the lattices of A and 
B. The cut in B is here along (110) and (110) directions. Pos- 
itive lobes of d-wave order parameter are shaded. Two qubits 
are shown, b) Version of a) using grain boundary (G) junc- 
tions. 

The dynamics of the device was considered in [Q. It 
is characterized by the phase difference ip between termi- 
nals A and B, which plays the role of the position of a 
quantum particle with mass M oc C, C being the clas- 
sical capacitance of the small terminal, in an effective 
two- well potential U(ip) (Fig. 2). It is the crucial advan- 
tage of clean DXD junctions, that the equilibrium phase 
iifo continuously depends on the angle between crystal 
lattices of A and B (and therefore on the d-wave order pa- 
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rameters in these terminals) in the interval [0, ir[ allowing 
for exponentially wide tuning of the tunneling rate [j7],|ll| . 
Moreover, due to time-reversal symmetry breaking in the 
system, states with tp = — tpa and tp = tpo are always de- 
generate and can be used as basic |0) and |1) states of a 
qubit 




FIG. 2. Effective potential profile of the system. Minima 
at ±ipo correspond to "up" and "down" pseudospin states of 
a qubit. The mismatch angle is fl = tt/8. 

The basic operations on a qubit are initialization, logi- 
cal operations (quantum gates) and measurement. Mea- 
surement is a two-step procedure and can be performed 
simultaneously on all qubits or selectively on individual 
or groups of qubits. The first step, collapse of the wave 
function, is achieved by grounding terminal A. Readout 
is facilitated by the existence of small persistent currents 
and magnetic fluxes (<C $o) which flow in opposite direc- 
tions in the |0) and |1) states 0,0. While too small to 
lead to unwanted inductive coupling between the qubits 
or decoherence, they can still be used to read out the 
state of the qubit once it was collapsed in one of the 
states with iipo, e.g. using a magnetic force microscope 
tip (which is removed during the computations). The es- 
timated magnetic moment of order 10 5 to 10 6 /i_e is on 
the resolution limit of commercial magnetic force micro- 
scopes. The same property can be used to initialize indi- 
vidual qubits or whole registers, since this small coupling 
to an external field can put the qubit in a desired (|0) or 
|1)) initial state. 

Let us now describe how logical operations can be re- 
alized in this system. In order to maintain coherence, the 
qubit 's electrodes A are isolated from ground while per- 
forming logical operations. The basic one-qubit logical 
operations are rotations around the x and z axes, X(9) 
and Z(4>): 

X{9) = e"^/ 2 ; (1) 
Z{4>) = e-™^' 2 . (2) 

Operation X(9), where 9 = 2tA/h and A are the tun- 
neling matrix elements, is provided by natural quantum 



beats between the two basis states |0) and |1). On the 
other hand, an effective rotation around the z axis is re- 
alized by lifting the degeneracy of the basis states by an 
amount exceeding the tunneling width. Thereby tunnel- 
ing between the basis states is suppressed and the nat- 
ural oscillations between the basis states, X(9), do not 
interfere with Z(<p) operations. The degeneracy between 
up/down states can be lifted in various ways. For exam- 
ple, it can be achieved by directly applying a localized 
magnetic field using a magnetic scanning tip. Other im- 
plementations will be discussed elsewhere. 

As stated, the idle-state of this system corresponds to 
the logical operation X(9). For a single qubit 'computer' 
this poses no problem as logical operations would be ap- 
plied sequentially without waiting times ("do-nothing" 
periods). In the case that a "do-nothing" period is de- 
sired, one can choose this time to be a multiple of the os- 
cillation period. Thus, using the above convention, this 
is equivalent to applying X(2nit) — 1, with n an integer 
and 1 the identity operator. The situation with more 
than one qubit is less straightforward. Here, we explic- 
itly need passive qubits (qubits which undergo no logical 
operations) to be "frozen" during operation on the active 
qubits (qubits over which a logical operation is applied). 
For instance, if Z((j>) is applied on qubit one, the state 
of passive qubits must not change during this operation. 
Since the application time of logical gates will typically 
be incommensurate with the time required for X{9) to 
be equal to 1, a scheme to freeze passive qubits is nec- 
essary. For this sake, it can be advantageous to have an 
idle-state where the energy of |0) and |1) are degenerate 
and tunneling is coherently blocked. One way to do this 
would be to temporarily enlarge the capacitance of elec- 
trode A by linking it with an external circuit as it was 
suggested by Ioffe et. al. in their "quiet-qubit" proposal 
However, such an approach brings the risk of losing 
coherence due to inelastic processes in the external nor- 
mal circuits. On the other hand, making the external ca- 
pacitor superconducting would bring uninvited evolution 
due to Josephson coupling between the external capacitor 
and electrode A. 

Our suggestion is to employ instead a technique of 
dynamic global refocusing closely related to refocusing 
methods of NMR 1 12 lj| and strong focusing of accelera- 
tors physics |0 . It relies on periodic perturbation of the 
two-well potential with amplitude 5E slightly exceeding 
the tunneling width. In this scheme, the energy shift be- 
tween the basis states is periodically varied from — SE to 
SE. Explicitly, this corresponds to the pulse sequence 



- Z{5Et/K) - Z{-8Et/K) - Z(6Et/U) ~ 



(3) 



This results in a time dependent angle of rotation 
around the z axes which is given, in the ideal case, by a 
triangular function of period 2r, the period of the refo- 
cusing sequence. 
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The evolution operator for a single qubit is then given, 
without approximation, by the Magnus expansion Jl3[ ] : 
U(t) = exp[-icr 2 J * dt'4>'(t')/2] = exY>[-i<7 z (SET-5ET + 
SEt — ■ ■ -)/2h] so that, in the worst case, it is equal to 
exp[±ia z $ET /2h]. For r sufficiently small this reduce 
to U(t) « 1. Hence, this yields a true idle-state as the 
information encoded by the qubits is not perturbed by 
tunneling nor by accumulation of relative phase between 
the basis states. The characteristic time scale of the refo- 
cusing pulse must be much less than the tunneling time 
(estimated in jjj as ~ 10 -8 s). 

It was recently demonstrated (Viola and Lloyd [|15| us- 
ing the spin-boson model; Viola, Knill and Lloyd [ [i6|Jl7| 
under more general assumptions) that in the limit of very 
small t, global refocusing leads to decoherence suppres- 
sion in the a x and a y channels (phase decoherence) pro- 
vided that cj) — 7r and that delays between the refocusing 
pulses are smaller, or of the order of, the correlation time 
of the environment [ft8f . This correlation time is given by 
the inverse of a natural cutoff frequency t c ~ cj~ and 
determines the fastest time scale of the environment. In 
the case of semiconductor-based structures, where deco- 
herence is due to phonons, t c is given by the inverse of 
the Debye frequency a;" 1 « 10~ 13 s p6[ . In the present 
situation, for r to be very small requires r <C tb, where 
tb <~ l/vf is the ballistic time (the time required for the 
formation of Andreev levels in the normal part of the 
system), I the size of the system and Vf the Fermi veloc- 
ity. Taking I ~ 10 3 A and Vf ~ 10 7 cm/s[[?|, we arrive 
at t < 10~ 12 s, a similar estimate as in |16 . Another 



potentially dangerous source of decoherence comes from 
the localized degrees of freedom (nuclear spins and para- 
magnetic impurities ) p9| ]. The estimates based on the 
central spin model |19[ show that the relevant energies 
correspond to much longer times, in excess of 10~ 8 s. 
(The same estimate can be made for the decoherence 
time from these subsystems.) On the other hand, the 
dynamics of a spin bath is much more complicated than 
the one of oscillator bath or spin boson models, and its 
behavior under global refocusing should be a subject of 
special investigation. 

Logical gates can be performed simultaneously with 
global refocusing pulses. Indeed, because the refocus- 
ing pulses obviously commute with Z '(</>), refocusing can 
be applied to all qubits (actives and passives) while per- 
forming Z((j)) on a qubit or in parallel on a group of 
qubits. The evolution of the active qubits is then given 
by exp[-ia z (J* dtp (?) + <!>)/ 2] « Z{4>). As a result, ap- 
plication of Z{<j)) on, e.g. the first qubit, in combination 
with the refocusing sequence yields the desired overall ac- 
tion on all qubits : Z(0)(g)l®- • -(g) 1. On the other hand, 
applying X(9) reduces to stopping refocusing pulses on 
the active qubits for a determined period of time. This 
also yields the desired overall action |20|. 

In order to create entangled states, non-local gates 



are required. Such an entangling two-qubit operation 
is realized in this system by opening the parity key join- 
ing two adjacent qubits, Fig 1. With this parity key 
open, a Josephson current flows between states of oppo- 
site phases. Thus, the combination |00) and |11) carries 
no current while |01) and 1 10) does. As a result, states 
of opposite phase will differ from those of identical phase 
by a Josephson energy Ej ~ 1 — cos(2<y9 )- The evolu- 
tion of a pair of qubits in this situation then corresponds 
to a conditional phase shift (CP) and, to an irrelevant 
phase factor, can be represented in the computational 
basis { |00), |01), 1 10) , 1 1 1) } as 



CP(j) = Diag(e I7/2 ,e-^/ 2 ,e-^/ 2 ,e^ /2 ), 



(4) 



with 7 = Ejt/h. Because CP(j) is diagonal in the com- 
putational basis, it commutes with Z(ifi). As a result, and 
under the assumption that the Josephson energy only 
weakly perturbs individual two- well potentials 21 1, the 
latter operation can be performed simultaneously with 
the global refocusing sequence. This condition can al- 
ways be realized by tuning the gate voltage on the parity 
key, thus varying its transparency and Josephson energy. 

Using the three basic operations defined above, it is 
possible to construct a Controlled-Not gate. This oper- 
ation, denoted CNij where i and j are the control and 
target qubits respectively, acts as : CN 12 \i, j) = \i, 
with © denoting addition modulo 2. Using the above ex- 
pressions for one- and two-qubits gate, CNi 2 is realized 
in this system, up to an irrelevant global phase factor, by 
the following sequence 

CN 12 = e l5 ^ 4 V 2 (7r /2)Z 2 (ix /2)A 2 (tt /2)Z 2 (ix /2) 

Z 1 (-k/2)CP(-k/2)X 2 (-k/2)Z 2 (-k/2)X 2 (-k/2). (5) 

In this expression, Xi(9) (Zi((j))) applies X(6) (Z((j>)) on 
the i th qubit while leaving the others unchanged (e.g., 
Zi(4>) = Z{4>) <g> 1 <g> • • • 1). 

In the setup of figure 1, it is possible to apply two-qubit 
gates only to adjacent qubits. It is therefore necessary 
to introduce a swap operator, denoted SWij , which ex- 
changes the states of qubits i and j. A swap on two 
adjacent qubits is realized by the following combination 
of Controlled-Not gates 



SW12 = CWi 2 CA 21 CA 12 . 



(6) 



Using this operator repeatedly, it is then possible to jux- 
tapose any chosen pairs of qubits and, as a result, to ap- 
ply Controlled-Not gates on any chosen pairs of qubits. 

Because of the commutation relations between the 
Pauli operators, combinations of rotations around the x 
and z axes generate SU(2), the group of 2 by 2 unitary 
matrices with determinant +1. Thus, it is possible to 
realize all one-qubit gates on this system. Furthermore, 
as been shown by Barenco et. ai, the set of all single 
qubit gates and the Controlled-Not is complete for quan- 



tum computation 22 1 . It is therefore possible to generate 
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all of SU(2 n ) with proper sequences of gates in such a 
n-qubit DXD superconducting quantum computer. 

In conclusion, we have shown that a solid state super- 
conducting quantum computer suggested in allows ap- 
plication of a complete set of quantum logical gates and is 
therefore a realization of a universal quantum computer. 
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